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ABSTRACT 

The recently developed universal form of the two-body 
problem involves several transcendental functions. Since 
these functions are evaluated so frequently, it is worthwhile 
to develop approximations that minimize the number of 
arithmetical operations required. This paper presents sev- 
eral such approximations, based on theories of Chebyshev 
and Knuth, with bounds for the errors incurred when using 
them. 
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OPTIMAL COMPUTING FORMS FOR THE 
TWO-BODY C AND S SERIES 


by 

C. R. Herron, E. R. Lancaster, andW. R. Trebilcock 
Goddard Space Flight Center 


INTRODUCTION 

The classical solutions of the two-body problem separate naturally into the three cases of el- 
liptic, parabolic, and hyperbolic motion, the mathematics being considerably different for each 
case. A unified formulation is possible, valid for all three cases, if certain transcendental func- 
tions, which we call the c and s functions, are introduced. 

The unified formulation is fully developed by Battin (Reference 1) and will not concern us here. 
The purpose of this paper is to present approximations for the C and s functions and their deriva- 
tive functions, which reduce significantly the computation times required for their evaluation when 
compared to those required by Taylor series expansions. 


THE C AND S FUNCTIONS 

The C and s functions are defined by 

S(x) = (x 1/2 ~ sinx 1/2 |/x 3/2 , x > 0 (1) 

= [sinh (-x) 1/2 - (~x) 1/2 ]y^(~x) 3/2 , x < 0 ; (2) 

C(x) = (l - cos x 1/2 )y/x , x > 0 (3) 

= [l - cosh (~x) 1/2 ]^x , x < 0 . (4) 


Since these functions are indeterminate for x = 0 and present accuracy problems when evalu- 
ated in the neighborhood of x = 0, it is natural to replace the above forms by the following series, 
convergent for all values of x : 


S(x) 


L 

i = n 


(-X ) 1 , 

( 2i + 3)! 


(5) 
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For large values of x, the convergence of these series will be slow. It is then convenient to 
use the following reduction formulas, easily derived from Equations 1 through 4: 


A(x) - 1 - xS(x) , 

( 7 ) 

2C(4x) = [A(x)] 2 , 

( 8 ) 

4S(4x) = S(x) + A(x) C(x) . 

( 9 ) 


THE C’ AND S’ FUNCTIONS 

The derivatives S' (x) and c' (x) are needed for certain problems of orbit determination, guid- 
ance, and optimization. From Equations 1 through 4 we obtain 


S' (x) - [C(x ) -3S(x)]y / 2x , 

C' (x) = [A(x) - 2C(x)]^2x . 

These forms suffer accuracy problems in the neighborhood of x = 0, again forcing us to series 
representations. Differentiating Equations 5 and 6, we have 


S' (x) 


L 

i = i 


( 2i + 3)! 


C' (x) 


□ 


E 

; = i 


( 2i + 2)! 


( 10 ) 


( 11 ) 


convergent for all values of x. 

For large values of x, the following reduction formulas (obtained by differentiating Equations 7, 
8, and 9) are useful. 


B(x) = 

S(x) + xS' (x) , 

(12) 

C* (4x) 

= -A(x)B(x) . 

(13) 

4S' (4x) = S' (x) 

+ A(x) C' (x) - B(x) C(x) . 

(14) 
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THE FIKE-KNUTH ALGORITHM 


Our first step in obtaining economical computing forms for Equations 5, 6, 10, and 11 was 
the construction of sixth degree polynomial approximations on various intervals in the sense of 
Chebyshev. In other words, these polynomials minimize the magnitude of the maximum error on 
the interval. The program to accomplish this was written by the third author, based on ideas of 
Stoer (Reference 2). The coefficients of these polynomials are given in the section entitled Numer- 
ical Results. 

Assume that the approximating polynomial has the form 

P(x) = a 0 +a 1 x + a 2 x 2 + a 3 x 3 + a 4 x 4 + a 5 x 5 +a 6 x 6 . (15) 

The evaluation of Equation 15 by the usual method of nested multiplication requires six multi- 
plications and six additions. However, by using recently developed polynomial evaluation methods 
(References 3 and 4), Equation 15 can be evaluated with four multiplications and seven additions. 
The form and parameters for the algorithm, as it applies to our functions, are given in the Numer- 
ical Results section. 

In the following description of the algorithm, a 6 is assumed to be positive. If a 6 is negative, 
a minor change is necessary. 

Fike f s modification of Knuth T s method begins with a conversion: let^ = y^, and let c k 3 a k ^. k 
for k = 0, 1, • • •, 5. Then compute 

P = 1 ( C 5 * 0 ' D " " C 2 " PC' • 

B' = C 4 - p(p + l) , E' = 2D' - B' + 1 , 

C' = c 3 - pB' , E" = 2D" - B' D' - C' , 

D' = p - B' , E'" = c, - B' D" . 

Find a real root q of the cubic equation* 

2q 3 + E' q 2 + E" q + E'" = 0 (16) 

•See Appendix A. 


3 


and compute 


A = \ B' - q , 

C = p - 2A , 

B = q - 2AC - A 2 , 

D = C' - q(l + D' ) - q 2 - D" - A 2 (1 + C) - BC , 

E = q 2 + qD' + D" - (A 2 +B)C , 

F - c Q - (q 2 + qD' + D" ) [c' - q( 1 + D' ) - q 2 - D"] . 

Then our polynomial can be evaluated as follows: 

r mx , 

q 2 = (q, + A) 2 , 

^3 = K + B ) K +C ) • 

p ( x ) = (q 2 + + D) (q 3 + E ) + F . 

In case a g < 0, let T(x) = -P(x) and perform all the steps above, except the last, for T(x). The last 
step should be 


P(x) - - [T(x)j - (q 2 + q 3 +D) (-q 3 -E ) - F . 

If the machine being used has a "load negative” feature that is equivalent in execution time to 
M load positive,” and if subtraction is likewise equivalent to addition, then this modification is equiv- 
alent to the original. 

As Fike points out, his method is a slight variation of that of Knuth (Reference 4), and since 
Knuth’s method was inspired by Motzkin (Reference 5), the three types bear a strong family re- 
semblance. Each begins with a polynomial in the form of Equation 15 with a 6 = 1, and solves for 
the parameters in the final evaluation scheme by expanding the scheme into a sixth degree poly- 
nomial and equating its coefficients with those of Equation 15. To admit treatment of the general 
polynomial of degree six, however, some transformation must be made so that a 6 =1. The most 
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straightforward way is making 

Q(x) = P(x)/a 6 

and then applying any of the three methods to Q(x), adding an extra step at the last in multiplying 
the result by a 6 . Fike specifies a different sort of transformation; his may be thought of as con- 
verting Equation 15 into 

a s a 4 a 3 a 2 a i 

R(x) = X + — X 5 + -7 X 4 + — X 3 + — x 2 + — X + a n . 

(X ' \jt fjr /jl z u 

Again, any of the three methods apply to R(x), and values of P(x) are obtained by using wx in the 
scheme for R(x), since R(^x) = P(x). 

This transformation, though a bit more complicated, is admirably suited to our particular prob- 
lem. The type of polynomial with which we are dealing has the not uncommon characteristic that 

I a 6 l < | a 5 | < • • • < | a 0 l 

and, in addition, |a 6 | is very small. For example, suppose the coefficients of form (15) are 

a 6 = + . 11 X 10‘ 10 . 

a 5 = - .21 x 10' 8 , 

a 4 = + .28 x 10' 6 , 

a 3 = - .25 x 1CT 4 , 

a 2 = + . 14 x 10" 2 , 

a j = - .42 x 1CT 1 , 

a Q “ + . 50 . 

If we use the division transformation, the coefficients b. of Q(x) are 

b 6 = + 1.0 . 

b 5 = - .18 x 10 3 , 
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+ .24 x 10 s , 


b 3 = - .22 xlO 7 , 

b 2 = + .12 xlO 9 , 

b x = - .36 x 10 10 , 

b 0 = + .44 x 10 11 . 


Here the errors in the numbers a. have become greatly magnified; worse yet, the arithmetic of 
parameter production using the large numbers is likely to suffer the effects of large error prop- 
agation. In contrast, Fike's transformation gives us 


C 5 = - .27 xlO 1 , 

c 4 = + .54 x 10 1 , 

c 3 = - .73 x 10 1 , 

c 2 = + .62 x 10 1 , 

c 1 = - .28 x 10 1 , 

c 0 = .50 . 


These numbers of manageable size lend themselves very well to whichever scheme we choose. For 
comparison, the two transformations above were evaluated by the Knuth algorithm for 40 points 
over the interval [- 1, +1], and the differences between these values and the true values of the 
polynomial were obtained. For the division transformation, the absolute value of the maximum 
error was .92 x 10~ 12 ; for the Fike transformation, this was .16 x 10 ~ 14 , a reduction by a factor 
of more than 500. Several other test cases were run, with results that apparently verify the conclu- 
sion that the Fike transformation used on this type of polynomial has a very definite advantage. 
There are, of course, other transformations that produce polynomials in which a 6 = 1. In general, 
one should use the transformation that keeps the coefficients of the transformed polynomial as small 
as possible. 
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NUMERICAL RESULTS 


The four approximation polynomials were generated for each of the intervals [- 1, +1], C- 2, + 2], 
[- 4, + 4], [- 16, + 16], converted to Equation 15 and parameters for the Fike evaluation scheme 
were obtained. In«each case, the values given by the final scheme were tested against "true" values 
of the original function for all multiples of .002 in the interval concerned. The true values came 
from expanding the power series of the function (1) for enough terms to guarantee that the relative 
error from truncation would be less than 10' 15 . The following tables exhibit, for each of the six- 
teen functions considered, the coefficients a £ for Equation 15, the parameters A, B, C, D, E, and F 
for the Fike scheme, and the maximum absolute errors for both methods. For comparison, a de- 
gree 4 approximation polynomial was evaluated by both methods for the functions C(x) and S(x) on 
the interval [- 1, +1], and the results are presented here also. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, March 27, 1968 
188-48-01-05-51 
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■mill 


C(X) = a 0 + a ? x + a 2 x 2 + a 3 x 3 + a 4 x 4 + a 5 x 5 + a 6 x 6 on jjh, hj 



a 0 = +0.4999999999999998 X 10° 
a, = -0.4166666666667176 X 10' 1 
a 2 = +o. 1388888888888999 X 10~ 2 
a 3 = -0.2480158725995993 X 10' 4 
a 4 = +0.2755731917059028 X 10' 6 
a 5 = -0.2087759200397967 X 10' 8 
a 6 = +0.1147134108311665 X 10' 10 

M a = 0.763 X 10' 15 


a Q =+0.4999999999999993 X 10° 
a 1 =-0.4166666666700118 X 10' 1 
a 2 = + 0.1388888888892785 X 10" 2 
a 3 = -0.248015866324180 7 X 10" 4 
a 4 =+0.2755731881710992 x 10" 6 
a 5 = -0.2088010277268315 X 10" 8 
a 6 =+0.1147215380312168 X 10" 10 

M a = 0.956 X 10" 13 


h = 4 

a 0 = +0.4999999999998401 X 10° 
a, = -0.4166666668808485 X 10' 1 
a 2 = +0. 1388888889138842 X 10' 2 
a 3 = -0.2480157659289839 X 10' 4 
a 4 = +0.2755731272513377 X 10' 6 
a 5 = -0.2089014196095935 X 10" 8 
a 6 = +0.1147636934013430 X 10' 10 

M a = 0.122 X 10' 10 


h = 16 

a 0 = +0.4999999894793170 X 10° 
a, = -0.4166675473500692 X 10' 1 
a 2 = +0. 1388889916034137 x 10' 2 
a 3 = -0.24798836331 19184 X 10' 4 
a 4 = +0.2755565077419917 x 10' 6 
a 5 = -0.2109148028487573 x 10" 8 
a 6 = +0. 1156091702389399 x 10" 10 

M a = 0.202 x 10~ 6 



q 2 = (9i + A) 2 

93 = (92 + B) ( 9 1 + C) 

C(X) = (q2 + q3 + D)(+q3 + E)+F on h, hj 
h = 1 h = 2 


A = +0.4513408582627891 X 10° 

B = +0.3744865190483202 X 10 1 
C = -0.2769272800754423 X 10 1 
D = + 0.9433565393074166 X 10 1 
E = + 0.1055413968372178 X 10 2 
F = +0.6288190624578802 X 10' 2 

A = +0.4512008438957284 X 10° 

B = + 0.3743936586512442 X lo 1 
C = -0.2769076432349482 XlO 1 
□ = + 0.9429681 327692907 X 1 0 1 
E = +0.1055053194443676 x 10 2 
F = +0.6284207351324604 X 10" 2 

M = 0.583 X10' 14 

M = 0.994 XlO" 13 

h = 4 

h = 16 

A = + 0.4507019590625572 x 10° 

B = +0.3740448131455307 x 10 1 
C = -0.2768317205414987 x 10 1 
D = +0.9414899439055362 X 10 1 
E = +0.1053701311149719 X 10 2 
F = +0.6272339359526764 X 10' 2 

A = + 0.4405766736959988 X 10° 
B = + 0.3669989101432331 X 10 
C = - 0.2752824996097087 x 10 1 
D = +0.9117484138879304 X lo' 

E = +0. 1026464040151929 X 10 2 
F = +0.6161613003116790 x 10' 2 

M= 0.122 XlO' 10 

M = 0.202 XlO' 6 
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on [-h, h] 


S(X) = a 0 + a| x + a 2 x 2 + a 3 x 3 + a 4 x 4 + asx 5 + a^x 6 


h = 1 

h = 2 

a 0 = + 0. 1666666666666665 X 10° 
a, = -0. 8333333333333568 X 10' 2 
a 2 = +0. 1984126984129264 X 10" 3 
a 3 = -0.2755731919939401 X 10“ 5 
a 4 = + 0.2505210785999854 X 10' 7 
a 5 = -0. 1605953765319026 x 10" 9 

a 6 = + 0. 7650283228592385 X 10“ 12 

a 0 = + 0. 1666666666666663 X 10° 
a ] = - 0. 8333333333352921 X 10” 2 
a 2 = + 0. 1984126984129057 X 10' 3 
a 3 = -0.2755731883077317 x 10" 5 
a 4 = + 0.2505210816170333 x 10" 7 
a 5 = -0.1606101133600677 x 10 " 9 
a 6 = +0.7647926042737674 X 10" 12 

M a = 0.555 X 10" 16 

M a = 0.566 X I0' 14 

II 

SZ 

h = 16 

a 0 = + 0. 1666666666666581 X 10° 
a, =-0.8333333334593103 X 10” 2 
a 2 = +0. 1984126984258407 X 10“ 3 
a 3 = -0.2755731292513204 x 10“ 5 
a 4 =+0.2505210496761327 x 10' 7 
a 5 = -0.1606691704048320 X 10' 9 
a 6 = + 0.7650122280184766 X 10" 12 

a 0 = + 0.1666666661133027 X10° 
a, = -0.8333338509758059 X 10" 2 
a 2 = + 0. 1984127524406629 x 10" 3 
a 3 = -0.2755570214946503 x 10" 5 
a 4 = + 0.2505123071826789 x 10" 7 
a 5 = -0.1618528418504030 x 10" 9 
a 6 = + 0.7694603615375217 x 10" 12 

M a = 0.720 X 10” 12 

M a =0.118 X IQ" 7 
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9 i = x 

q 2 = (q, + a ) 2 
^3 = (^2 + b) (q, + c) 

S(X) = (q 2 + q3 + D) ( + q 3 + E) + F on [- h, hj 
h = 1 h = 2 


A = + 0. 1030541 H 0544949 X 10° 
B = +0. 1357446199107850 X 10 1 
C = -0. 1709888012144409 X 10 1 
D = +0. 1803960616654583 X 10 1 
E =+0.2062171852872241 x 10 1 
F = +0.2130010466488949 x 10" 1 

M = 0.236 xiO" 14 


A = + 0. 1028274494783523 X 10° 

B = + 0. 1356879505417743 XIO 1 
C = - 0. 1709784631095070 X 10 1 
D = + 0. 1802832347589989 X 10 1 
E = + 0.2060949727430426 X 10 1 
F = + 0 . 2 1 28753997322974 XIO -1 

M = 0.808 XIO " 14 


h = 4 

A = + 0.1026453527945748 x 10° 

B = + 0.1356011711587295 x 10 1 
C = -0.1709549340846876 x 10 1 
D = + 0.1800557277079013 x io 1 
E = +0.2059246874696585 x 10 1 
F = + 0.2125209137884825 x 10" 1 

M = 0.722 XIO -12 


h = 16 

A = + 0. 9898746283297480 XIO" 1 
B = + 0. 1338586121335949 X 10 1 
C = -0. 1704756186376375 X 10 1 
D = + 0. 1754906650979839 x 10 1 
E = + 0.2025050653157090 x 10 1 
F = + 0. 2056593593968585 x 10" 1 

M= 0.118X10" 7 
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II I II ■■ IB I 


■ III I I II 


C'( X) = a 0 + a] x 


+ a 2 x 2 + a 3 x 3 + a 4 x 4 + a 5 x 5 + a 6 x 6 on 


[-h, h] 


h = 1 

a 0 = - 0.4166666666666430 X 10" 1 
a, = + 0.2777777778467318 X 10' 2 
a 2 = -0. 1488095238678453 x 10" 3 
a 3 = + 0.6613751098214413 x 10" 5 
a 4 =-0.2505208412532178 x 10' 6 
a 5 = + 0.8269965205281566 x 10" 8 
a 6 =-0.2412214891589441 x 10" 9 

M a = 0.217 X 10 -16 


h = 4 

a 0 = -0.4166666651 125364 X 10' 1 
a, = +0.2777780643726382 X 10' 2 
a 2 = -0. 1488097663598733 X 10" 3 
a 3 = +0.6612326049924104 X 10' 5 
a 4 = - 0.2504581368842600 X 10' 6 
a 5 = + 0.8437138609417991 X 10' 8 
a 6 = -0.2463133603265637 X 10" 9 

M a = 0.320 x 10' 12 


h = 2 

a 0 =-0.4166666666606742 X 10" 1 
a 1 =+0. 2777777822034584 X 10' 2 
a 2 =- 0. 1488095275535492 x 10" 3 
a 3 =+0.6613668137463213 X 10" 5 
a =-0.2505171918626241 X 10" 6 
a = + 0 . 8303 1 36595603477 x 1 0 " 8 
a 6 =-0.24223 1668 1583509 X ]0" 9 

M a = 0.251 X 10" 14 


h = 16 

a 0 = - 0.4166666641763858 X 10" 1 
a, = + 0.2777780078584500 X 10" 2 
a 2 = -0.7440478621862503 X 10" 4 
a 3 = +0. 1102220894089232 X 10" 5 
a 4 = - 0. 1043798352532477 X 10" 7 
a 5 = + 0.6938557071056230 X 10"'° 
a 6 =-0.3366982823031963 X 10' 12 

M a = 0.527 X 10~ 8 
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I 


9l = v^96 X 

q 2 = (^i + A ) 2 

93 = (q 2 + b) (q i + c) 

C'( X) = (q2 + q 3 + D) ( - q3 - E) - F on [- h, hj 


h = 1 

A = + 0.642 1 500675794880 x 1 0 " 1 
B = +0.9631412054424315 x 10° 
C = -0.1485378935681960 X 10 1 
D = + 0. 1206516480446427 X 10 1 
E = + 0. 1262370275488431 x 10 1 
F = + 0.2235785774036898 x 10 " 2 

M = 0.179 x 10" 15 


h =4 

A = +0.6405925631551870 x 10 _1 
B = +0.9624584033933687 X 10° 
C = -0.1485269967153406 x 10 1 
D = +0.1205018912730453 x 10 1 
E = + 0.12613947171 19364 X 10 1 
F = + 0.2210887038704210 x 10" 2 

M = 0.320 x IQ" 12 


h = 2 

A = + 0.6415596248585610 X 10 -1 
B = + 0. 9629599640982438 X 10° 
C = - 0. 1485350658374242 X 10 1 
D = +0. 1206153148806730 X 10 1 
E = +0. 1262088821410514 X 10 1 
F = +0.2230745797490601 x 10 _2 

M = 0.268 X IQ" 14 


h = 16 

A = + 0.6208183431484613 X 10 
B = + 0.9523193847619871 X 10° 

C = -0.1483582934430130 X 10 1 
D = + 0.1 182131 124018069 X 10 1 
E = + 0. 1247277466997771 X10 1 
F = + 0. 1829570481293727 X 10 ' 2 

M = 0.527 x 10" 8 
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S'(X) = ag + a, x + a 2 X 2 + 


on jj- h, hj 


a 3 x 3 + a 4 x 4 + 85 x5 + x 6 


h = 1 

a 0 = -0.8333333333333210 x 10 “ 2 

a, =+0.3968253968616768 x 10 ~ 3 

a 2 = - 0. 1653439153716376 X 10 “ 4 

a 3 = +0.6012503110063301 x 10 “ 6 

a 4 = -0.1927084107177350 x 10 “ 7 

a 5 = + 0.5511761621292874 x 10 " 9 

a = -0.1418571972378231 x 10’ 10 
6 

M a = 0.260 X 10 “ 17 


h = 4 

a 0 = -0.8333333325953303 x 10 “ 2 
a, = +0.3968255472561213 x 10 ’ 3 
a 2 = -0. 1653440305444520 x 10 “ 4 
a 3 = +0.601 1754909568150 x 10 “ 6 
a 4 = -0.1926786201990815 x 10 “ 7 
a 5 = +0.5599576811597955 x 10 “ 9 
a 6 = -0.1442776137237285 X 10" 10 
M a =0.168 x 10 “ 13 


h = 2 

a 0 = -0.8333333333304809 x 10 ~ 2 
a, = +0.3968253991531185 x 10" 3 
a 2 = -0. 1653439171256403 x 10“ 4 

a 3 = + 0.6012459474285134 x ] 0~ 6 
a 4 = -0. 1927066737597982 x 10“ 7 
a 5 = +0.5529210296030711 x 10 " 9 
a 6 = -0.1423381311296310 x 10 "'° 

M a = 0.134 X 10 " 15 


h = 16 

a 0 = -0.8333333321480760 x 1CT 2 
a, = + 0.3968255178485000 X 10 " 3 
a 2 = -0.8267196924462379 X 10 “ 5 
a 3 = + 0 . 1 002046526767437 X 1 0“ 6 
a 4 = -0.8029333915166488 x 10“ 9 
a 5 = +0.4617817733427159 X 10“ 11 
a 6 = -0.1978183008506138 X 10“ 13 

M a = 0.277 X 10 “ 9 
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9 i = x 

92 = (9i + A) 2 

93 = (92 + B ) (9i + C) 

S'(X) = (92 + ^3 + D ) (“93 " E) - F on £- h, hj 


h = I 

A =-0.4272502910304863 x |0‘ 1 
B =+0.4460409572215307 X |0° 
C =-0.1020276052465536 X I0 1 
D =+0.4064530281034859 x 10° 

E =+0.3450015283629582 x 10° 

F =+0.2885054406289098 x |0“ 2 

M= 0.529 x 10" 16 


h =4 

A =-0.4274056995233864 x 10“ 1 
B =+0.4458621272744412 x 10° 

C =-0.1020307887595762 X I0 1 
D =+0.4061739766937458 x 10° 

E =+0.3449170623874602 X 10° 

F =+0.2876661053629565 x I0‘ 2 

M = 0.169 x 10” 13 


h = 2 

A =-0.4272336986330691 X|0 -1 
B =+0.4460090638507059 x |0° 
C =-0.1020282840851638 X I0 1 
D =+0.4064007086884238 X |0° 

E =+0.3449898120978949 X |0° 

F =+0.2883373149053849 x I0 -2 

M= 0.182 X 10 -15 


h = 16 

A = -0 . 43 1 0245077956 152 X 1 0" 1 
B =+0.4428994350437274 x 10° 
C =-0.1020789569928175 x 10 1 
D =+0.4016042135773923 x 10° 

E =+0.3434241762739658 X 10° 

F =+0.2744481265299470 x I0" 2 

M = 0.277 x 10" 9 
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P(X) = a 0 + ajx+a 2 x2 + agx3+a 4 x 4 


C(X)#4 [-l, + l] 

a 0 = + 0.5000000000007167 X 10° 
a } =-0.4166666601424471 X10" 1 
a 2 = + 0.1388888879568642 x 10" 2 
a 3 = -0.2480419696201834 x 10" 4 
a 4 = + 0.2755932664579228 X 10 -6 

M a = 0.130 X 1(T 9 


S(X)#4 [-1, + l] 

a 0 = + 0.1666666666667142 X 10° 
a, = -0.8333333283147393 x 10" 2 
a 2 = + 0 . 1984126977913694 X 10 -3 
a 3 = -0.2755932664326337 X ]0' 5 
a 4 = + 0.2505344666229896 x 10" 7 

M a = 0.100X IQ -10 


q i = Ax , 

q 2 = ( c h + B ) 2 

p (x) = ( q 1 + q 2 + C) ( q 2 + D) + E 


C(X) H 



S(X)^4 £l, +l] 


A = + 0.2291221893900772 x 10" 1 
B = -0.7655416156428518 x 10° 
C = + 0.2592589041826887 X10° 
D = + 0.4011554686880556 X10° 

E = -0.3345008393894142 X 10° 

M = 0.130X10 -9 


A = + 0.12581 04947765253 X 1 O’ 1 
B = -0.5959855810891701 x 10° 
C = + 0.1 104585254341674 X 10° 
D = + 0.2038526154314646 X 100 
E = - 0 . 9365973340739070 X 1 O' 1 

M = O.IOOXIO" 10 
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and since v and t have the same sign, we have 


V + 



> | 2 q+1 1 1 


This simply says that the length of the initial bracketing interval is less than, or equal to, the mag- 
nitude of the small end; in turn, this means that the large end of the interval is at most twice the 
magnitude of the small end. 

Now, consider how the endpoint values would be represented in floating-point binary arithmetic 
(normalized) with an r-bit fraction. If the difference between their binary exponents is at most 1 
(which is what we are getting at above), then it can be seen that the number of distinct points in the 
initial bracketing interval is at most 2 r . Therefore, the number of interval -halving iterations 
needed— that is, the number of times one reduces his choice of points in the interval by one-half— 
is at most r. Moreover, it often turns out that f is nearly (or exactly) zero at an end point of one 
of the half-intervals, so that r iterations are not always needed. 

We have treated the special case 1 1 1 ^ | v| , but we need not restrict ourselves to it. The number 
of interval-halving iterations needed depends upon the size of t and, if one is willing to iterate a 
bit more, he can find the initial bracketing interval more quickly by increasing t; the converse of 
this also holds. 
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